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\Q ■ Abstract. We investigate the pulsation-phase dependent properties of the circumstellar dust shell (CDS) of the 

OH/IR star OH 104.9+2.4 based on radiative transfer modeling (RTM) using the code DUSFY. Our previous 
study concerning simultaneous modeling of the spectral energy distribution (SED) and near-infrared (NIR) visi- 
bilities (Riechers et al. I2004t has now been extended by means of a more detailed analysis of the pulsation-phase 
dependence of the model parameters of OH 104.9+2.4. In order to investigate the temporal variation in the spatial 
structure of the CDS, additional NIR speckle interferometric observations in the K' band were carried out with 
• the 6 m telescope of the Special Astrophysical Observatory (SAO). At a wavelength of A = 2.12 /xm the diffraction- 

Oh' limited resolution of 74 mas was attained. Several key parameters of our previous best-fitting model had to be 

adjusted in order to be consistent with the newly extended amount of observational data. It was found that a 
5_l ■ simple rescaling of the bolometric flux Fboi is not sufficient to take the variability of the source into account, as 

| the change in optical depth r over a full pulsation cycle is rather high. On the other hand, the impact of a change 

in effective temperature T e ff on SED and visibility is rather small. However, observations, as well as models for 
other AGB stars, show the necessity of including a variation of T e s with pulsation phase in the radiative transfer 
models. Fherefore, our new best-fitting model accounts for these changes. 
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1. Introduction 

OH 104.9+2.4 (IRAS 22177+5936, AFGL2885, NSV 
25875, IRC0 243) is an OH/IR type II-A class star. Such 
objects are highly evolved asymptotic giant branch (AGB) 
stars exhibiting dusty envelopes, which are believed to be 
the main source of dust grains in our Galaxy. While the 
peak of the SED is expected to be around 6 - 10 /im for 
OH/IR stars, the SiO and Si02 features at 9.7/iim and 
18 /im are usually found in absorption, indicating a high 
optical depth of the CDS. This is consistent with the 
high mass- loss rates (10~ 7 - 10" 4 M Q /yr, Habing 1996) 
that have been measured for these long-period variables 
(LP V) . A large number of observational constraints were 
found for this class of stars, leading to detailed models of 
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the evolution of these objects and the properties of their 
CDS (see Riechers et al. 120041 hereafter Paper I, and 
references therein). For OH 104.9+2.4, an outflow velocity 
of v e — 15km/s was measured by te Lintel Hekkert et 
al. I[1991[) . and the absolute distance was determined 
to be D = 2.38 ± 0.24kpc (Herman & Habing [TM^i . 
From IR observations, we derived a pulsation period of 
P = 1500 ± lid, as well as a bolometric flux ranging 
from Fboi = 0.7 x 10 _10 W/m 2 at minimum phase to 
-Fboi = 2.3 x 10 -10 W/m 2 at maximum phase (see Paper 
I for details). 

The purpose of the present paper is to discuss 
our improved best-fitting radiative transfer model of 
OH 104.9+2.4. The paper is organized as follows: in 
Sect. 2, results of our new NIR speckle interferometric 
observations are presented. Sections 3.1 and 3.2 summa- 
rize the results of our best-fitting models from Paper I 
(hereafter model M-I) and our new best-fitting model 
(hereafter model M-II). Section 3.3 tests the model 
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M-II quantitatively by comparing it with the SED data 
collected in Paper I for different pulsation phases. Section 
4 summarizes and evaluates the results and closes the 
discussion with an outlook on future research. 

2. Speckle interferometric observations and data 
reduction 

In addition to our speckle interferometric observations in 
2002, which were presented in Paper I, K'-hand speckle 
interferograms of OH 104.9+2.4 were obtained with the 
Russian 6 m telescope of the Special Astrophysical 
Observatory on October 11, 2003. The data were recorded 
with our HAWAII speckle camera through an interfer- 
ence filter with a center wavelength of 2.12 //m and band- 
width of 0.21 (im (K' band). Additional speckle interfero- 
grams were taken for the unresolved reference star 2MASS 
J22191350+5949584. With a pixel size of 28.7 mas and 
seeing of 2.6arcsec, 4100 object frames and 4600 frames 
of the reference star were taken, each with an exposure 
time of 218 ms. These interferograms were used to com- 
pensate for the speckle interferometry transfer function. 
The iT-band visibility function of OH 104.9+2.4 was de- 
rived from the speckle interferograms using the speckle in- 
terferometry method (Xabevrie ll970|) . The reconstructed 
two-dimensional visibility at 2.12 /im and azimuthally av- 
eraged visibility profile are shown in Fig. 1. Given the ac- 
curacy and resolution of our visibility measurements, we 
find that the diameter ratio between major and minor axes 
differs from unity by less than 6%, which justifies the as- 
sumption of spherical symmetry. 

From the pulsation period derived for OH 104.9+2.4 
in Paper I, the phase of observations presented in this pa- 
per is determined to be <I> 2003 = 0.25, as compared to the 
phase of the 1996 ISO spectrum ($ ISO = 0.5) and to the 
phase corresponding to the 2002 SAO visibility measure- 
ments ($ 2002 = 0). Taking into account the uncertainties 
in the pulsation period caused by the low number of photo- 
metric data points and the resulting limitation to a rather 
simple fitting procedure, we estimate the total error of 
these phases to be ~0.02. 

3. Dust shell models 

3.1. Previous results: Model M-l 

As reported in Paper I, for the previous analysis of the ra- 
diative transfer in the CDS of OH 104.9+2.4, - 10 6 models 
were calculated. To provide a reference for our SED model 
computations we used the 1996 ISO spectrum, while infor- 
mation about the spatial structure was given by our 2002 
SAO if '-band visibility. To simultaneously fit the SED 
and visibility data (i.e. with otherwise identical model pa- 
rameters), -Fboi was adjusted for the epoch of the visibility 
measurements with respect to the SED model, since the 
ISO spectrum and the 2002 visibility data were taken at 
different pulsation phases of OH 104.9+2.4. Therefore, M- 
I consists of an SED model where the bolometric flux was 



calculated from the 1996 ISO spectrum and from a visibil- 
ity model using a bolometric flux for the 2002 SAO phase, 
which was derived from photometric measurements by dif- 
ferent authors (see Paper I, Sect. 3). With these observa- 
tional constraints and phase dependencies, we have found 
a model that fits the observations. The main parameters 
of M-I are summarized in Table ^ 

3.2. The new best-fitting model: M-II 

Based on the new constraints given by the 2003 SAO vis- 
ibility, some changes to M-I had to be introduced. The 
most significant change is that in addition to adjusting 
-Fboi, a phase-dependent optical depth had to be intro- 
duced in order to be able to explain the observations. This 
approach was recently suggested based on analysis of the 
depth of the 9.7 /itm feature for other OH/IR stars (see 
Suh 200H and references therein). In addition, the new 
model considers variations in the effective temperature of 
the central star, which is a less sensitive model parameter 
than the optical depth. All results for M-II are summa- 
rized in Table E| The final model is shown in Fig.[21 

3.2.1. Step 1: improved SED fit 

Apart from the goal of obtaining a good fit for multiple 
pulsation phases, another aim of the ongoing research was 
to find a dust type that is capable of reproducing, in par- 
ticular, the 18 fim absorption feature better than the dust 
composition used in M-I. Therefore, we included in our 
analysis several dust types not inherent to the DUSTY 
code (Ivezic et al. I1999|l in order to widen our parameter 
space grid. The only dust type that was found to repro- 
duce the SED of OH 104.9+2.4 significantly better than 
the dust composition used in M-I is the cold silicates from 
Suh l|1999j) . With this dust type, we have found a model 
with an Ivezic et al. (1999) radiatively driven wind that 
gives a slightly better SED fit than corresponding models 
with p(r) oc r~ 2 . In Fig.0 the SED in the 10 /im regime 
for M-I and the final M-II are shown for comparison. 

The figure clearly illustrates how this new model is an 
improvement in this region. At this point, we would like to 
note that the ISO SED of OH 104.9+2.4 shows some rather 
small emission features at 33.3 /im and 40.6 /im, which are 
not considered in detail by our model even with the up- 
dated (as compared to M-I) dust composition. Suh ( 2002 ) 
ascribes these features to crystalline silicates. However, it 
is difficult to extract information from those rather weak 
features and especially to positively identify them with 
specific dust components. Therefore, proper identification 
of these features may not be possible with current mod- 
els. But, since these features are fairly weak, accounting 
for them in more detail would not significantly alter the 
conslusions of our study. 

As we decided to use a model incorporating a ra- 
diatively driven wind, we had to comprehensively check 
our parameter sets to match the observed outflow veloc- 
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Fig. 1. Observations of 
OH 104.9+2.4 on October 
11, 2003, with the SAO 
6 m telescope. Left: 2- 
dimensional visibility 
function of OH 104.9+2.4 
at 2.12/xm. Right: 
Azimuthally averaged 
visibility of OH 104.9+2.4 
in the K' band. 



Table 1. The derived and adopted physical parameters of OH 104.9+2.4 as provided by Model M-I. 



Parameter 



Value 



Effective temperature (black-body) 
Temperature at inner CDS boundary 
Density profile within the CDS 
Relative CDS thickness 
Dust-grain distribution function 
Minimum grain size 
Maximum grain size 
Dust types 

Optical depth 



Radius 

Radius of inner CDS boundary 
Mass-loss rate 
Bolometric flux 



Tea = 2500 + 500 K 

r in = 1000 + 200 K 

p{r) <x r~" with n = 2.0 ± 0.1 

^ = 10 p with p = 5tT 

MRN, n(a) oc a~ 35 (Mathis, Rumpl, & Nordsieck[T577Jl 
a min = 0.005 ± 0.003 Mm 
a max = 0.2 ± 0.02 pm 

95 % warm silicates from Ossenkopf, Henning & Mathis ( 1992 ) 
5 % astronomical silicates from Draine & Lee 1119841 

T"0.55fim = 158 ± 7 
T2.2 M m = 6.5 ± 0.3 



Tg.TL 



14.0 + 0.1 



R+ ~ 600 R Q = 2.79 AU 
R in = 9.1 R* = 25.4 AU 
M = 2.18 x 10" 5 M Q /yr 
F bol = 0.7 x 10" 10 W/m 2 



ities. The new best-fitting model results in an average 
terminal outflow velocity of v^f~ u = 16.8 km/s, which is 
in fair agreement with the observed outflow velocity of 
v c — 15 km/s from te Lintel Hekkert et al. (1991). 

In addition, the effective temperature was adjusted to 
T e g = 2250 K for phase $ = 0.5. Considering the error 
bars found for M-I, this is only a minor adjustment. 

3.2.2. Step 2: fitting multiple pulsation phases 

Comparing the 2002 and 2003 SAO visibility data, it is 
obvious that a simple rescaling of -Fboi is insufficient for 
reproducing both observational constraints (SED and vis- 
ibility) with a single model. The visibility scales with 
V / 7tioi, and the difference in Fboi between the 2002 and 
the 2003 measurements is AF bo i = 0.8 x 10" 10 W/m 2 
(based on the assumption of a cosine-like pulsation phase, 
see Paper I), corresponding to 35% of . A differ- 

ence of this order of magnitude cannot be derived from 
the SAO observations by merely rescaling the bolometric 
flux. Therefore, at least one additional parameter has to 



be changed to account for the observed variation of the 
visibility. To explain the observed change of depth for the 
9.7pm absorption feature of other OH/IR stars, different 
authors have suggested a change in optical depth with 
pulsation phase (see Suh 120041 and references therein) . 
With both proper scaling of Fboi as a function of the 
pulsation phase (4> 2002 = 0, F^ 2 = 2.3 x 10- 10 W/m 2 ; 
$2003 = o.25 5 ir2003 = ]_ 5 x 1Q -io w/m 2 ) and a linear ad- 
justment of the optical depth (r^ = 8.5, rfi, ^ = 3.5, 
Tf.^m = 6.0), the 1996 ISO SED, as well as the 2002 
and 2003 SAO visibilities, can be reproduced by a single 
model, as shown in Fig.El 

Most detailed pulsation models for AGB stars predict 
an increase of T c r between minimum and maximum pulsa- 
tion phases. For Mira stars, Bessell et al. ( 1996J) predict the 
maximum in T c g to be at $ ~ 0.8. While radial pulsations 
of OH/IR stars are driven by the same mechanism, it is 
difficult to assess whether the same $-T e fj-relation is true 
for these objects. However, our primary interest is to find a 
suitable model that properly describes the relative changes 
between $ = 0.5, $ = 0, and $ = 0.25. Based on Bessell et 
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X [^ m ] spatial frequency [cycles/arcsec] 

Fig. 2. SED (left) and 2.12 /im visibility (right) of OH 104.9+2.4 for our best-fitting model M-II. In order to take 
into account the different epochs of the ISO (1996) and SAO (2002 & 2003) observations at 2.12 /im, bolometric flux 
Fboh optical depth r, and effective temperature T e e had to be adjusted. The effective temperature of the central star 
is T cff = 2250 K for $ = 0.5, T cS = 3150 K for $ = 0, and T e $ = 2800 K for $ = 0.25. F_i denotes the bolometric 
flux in 10 -10 W/m 2 , while r 2 .2 is the optical depth at 2.2 /im. 
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Fig. 3. SED of OH 104.9+2.4 in the mid-IR for models M- 
I (Paper I) and M-II (this work). Improvement in fitting 
the silicate absorption features is mainly due to the newly 
included dust composition. 

al. H199(i|l and taking into account that OH 104.9+2.4 is a 
cool star, we find T off = 2250 K for $ = 0.5, T cff = 3150 K 
for $ = 0, and T cS = 2800 K for $ = 0.25 are appropriate 
to account for the phase-dependent change of the central 
star's effective temperature. Nevertheless, a change in T e g 
has only a small impact on SED and visibility and almost 
the same effect as would a small adjustment of r in the 
opposite direction. 

As can be seen from Table [3] the maximum grain size 
imax was increased for this model as compared to M-I. 
The new value of a max = 0.28 /im for M-II deviates from 
flmax = 0.20 /im for M-I by 4 cr, but the change was nec- 
essary in order to enable us to scale the optical depth 
properly. This does not come as a surprise, since, chang- 



ing two different parameters in opposite directions may 
lead to nearly identical models, as pointed out in Paper I. 
Therefore, we were able to find different models of simi- 
larly good quality in Paper I. The new best-fitting model 
M-II is much more constrained by observations than M- 
I, as the third observational constraint turned out to be 
essential in order to show the necessity of adjusting the 
optical depth for different phases of OH 104.9+2.4. 

Concerning the 2003 SAO visibility, M-II cannot repro- 
duce the observations beyond 12 cycles/arcsec very well. 
This may be due to the large observational errors close to 
the diffraction limit, which are possibly slightly underesti- 
mated by the error bars shown here, as there is no obvious 
physical reason why the visibility curvature should change 
as drastically as suggested by Fig.|3in that regime between 
the dates of the two SAO observations. 

3.2.3. Further model results 

The top left panel of Fig.0] shows the fractional flux con- 
tributions of the emerging stellar radiation, scattered ra- 
diation, and thermal dust emission as a function of wave- 
length for model M-II. At 2.2 /im, the flux is dominated 
by scattered light for all phases (see Table |3J). The contri- 
bution of direct stellar light (attenuated flux) rises from 
^6% at minimum phase to ^25% at maximum phase. 

The top right panel of Fig. 01 shows the normalized in- 
tensity distribution at 2.13 /im as a function of angular 
distance. The barely resolved central peak corresponds to 
the central star. Here, the increase in the flux contribution 
of the central source towards maximum phase is reflected 
by the rise of the central peak, as it does indeed become 
larger. In addition, the inner rim of the CDS itself moves 
outwards. 

The bottom left panel of Fig.0] shows a plot of the 
dust temperature as a function of angular distance. The 
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Table 2. The derived and adopted physical parameters of OH 104.9+2.4 as provided by Model M-II. Parameters that 
were treated as not time-dependent are indicated with *. A distance of D = 2.38 + 0.24 kpc (Herman & Habing 1985) 
and an outflow velocity of v c — 15kms _1 (te Lintel Hekkert et al. 1991) were assumed. 



Parameter 






$ = 0.5 
(1996) 


<E> = 
(2002) 


$ = 0.25 
(2003) 


Black-body effective temperature 


T e ff 


[K] 


2250 


3150 


2800 


Temperature at inner CDS boundary* 


Tin 


[K] 


1000 






Relative CDS thickness* 


''out 
Hn 




10 5 






Minimum grain size* 




H 


0.005 






Maximum grain size* 


&max 


H 


0.28 






Dust-to-gas ratio* 






0.005 






Dust grain bulk density* 


r s 


[gem -3 ] 


3.0 






Optical depth 


TO. 55 |im 




146.7 


60.4 


103.6 




T2.2 |im 




8.5 


3.5 


6.0 








13.9 


5.7 


9.8 


Radius 


R+ 


[Re] 


729 


675 


691 






[AU] 


3.39 


3.14 


3.21 






[mas] 


1.43 


1.32 


1.35 


Radius of inner CDS boundary 


Rin 


[R*] 


8.3 


17.5 


13.7 






[AU] 


28.3 


55.0 


44.0 






[mas] 


11.9 


23.1 


18.6 


Mass-loss rate 


M 


[10- 5 M Q /yr] 
[10- 10 W/m 2 ] 
[10 4 L Q ] 


3.09 


5.69 


5.17 


Bolometric flux 


F hol 


0.7 


2.3 


1.5 


Luminosity 


L 


1.23 


4.03 


2.63 



Density profile within the CDS* radiatively driven wind (Ivezic et al. ^999 1 

Dust-grain distribution function* MRN, n(a) oc a -3,5 (Mathis, Rumpl, & Nordsieck I1977t 

Dust type* S99 cold silicates (SuhUl^ 



inner CDS radius n n can be derived directly from the ra- 
dial profile, as it is the point at which the dust reaches 
Ti n = 1000 K. It corresponds to an angular radius of 
■d in = H.9mas for $ = 0.5 (for M-I: tf in = 10.5mas), 
&- m = 18.6 mas for $ = 0.25, and i? in = 23.1 mas for $ = 0. 
This again shows that the inner rim of the CDS moves out- 
wards as the pulsation phase rises from minimum to max- 
imum, mainly a consequence of increasing L+. According 
to our model, the change in the inner dust shell radius be- 
tween minimum and maximum phases is approximately a 
factor of two. 

The bottom right panel of Fig.0] shows the total opti- 
cal depth as a function of wavelength. As can be expected 
for the chosen chemical composition, local maxima of t\ 
are found for the silicate absorption features at 9.7/im 
(SiO stretching vibrations) and 18 fim (Si0 2 bending vi- 
brations). The shape of the profile remains the same for 
all pulsation phases, but the optical depth is higher for all 
wavelengths towards the minimum pulsation phase. 

The stellar radius drops by 8% between minimum and 
maximum phases (from 729 to 675 i?©), which is con- 
sistent with the Stefan-Boltzmann law. The radius at 
$ = 0.5 is approximately 18% larger as compared to the 
value derived for M-I. As described in Paper I, the mass- 
loss rate M can be calculated using the observed outflow 
velocity v c = 15km/s and our derived luminosities for 
the different phases (see Table For minimum phase, 
M = 3.09 x 10 -5 M /yr is obtained, which is roughly 50% 
larger than the value found for M-I. For maximum phase, 



Table 3. The contributions to the total flux at 2.2 /im 
for different pulsation phases $ as given by Model M-II. 



component 


$ = 0.5 


$ = 


$ = 0.25 


attenuated 


5.9% 


25.1% 


15.0% 


scattering 


55.9% 


42.2% 


51.7% 


thermal 


38.2% 


32.7% 


33.3% 



we derive M = 5.69 x 10 5 M©/yr, which is by a factor of 
1.8 larger than the mass- loss rate at minimum pulsation 
phase. Heske et al. lfTMi|l found M = 5.58 x 10~ 5 M©/yr, 
which is in good agreement with our result. 

3.3. Quality of the SED modeling at different phases 

In Sect. 3.2 it was shown that model M-II is able to 
properly reproduce the ISO SED and the SAO visibili- 
ties. But, as argued in Paper I, the SED varies in over- 
all flux as well as in shape, depending on the pulsation 
phase. Therefore, this has to be tested if M-II is to re- 
produce the SED of OH 104.9+2.4 at different pulsation 
phases. Unfortunately, only few photometry is available 
besides the ISO spectrum at minimum phase; while for 
some data, the pulsation phase is uncertain, as the ob- 
serving date is not specified in the corresponding publica- 
tions. Therefore, we decided to perform only a qualitative 
test on the data used in Paper I to derive the SED at 
maximum phase. Figure [5] shows the results. The models 
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Fig. 4. Modeling results of OH 104.9+2.4 for our best-fitting model M-II at the pulsation phases of the observations 
(see Table EJ. Top left: Fractional contributions to the total flux. The contribution from direct stellar light is small 
compared to the contribution of dust scattering and thermal emission to the physical processes in the CDS (see Table 
03). Top right: The normalized intensity profile at 2.13/xm. The sharp central peak corresponds to the central source 
of radiation. Bottom left: Dust temperature as a function of angular distance. The point where Tdust = 1000 K 
indicates the inner radius of the dust shell. As the figure shows, due to the variability of the central source, the inner 
dust shell boundary moves from 11.9 mas to 23.1 mas from minimum to maximum pulsation phase. Bottom right: 
Wavelength depencence of the total optical depth. 



at minimum and maximum phase properly constrain the 
range of observational data. Thus, Model M-II seems to 
correctly reproduce the change in shape and the scaling 
with Fboi of the SED of OH 104.9+2.4. 

4. Summary and conclusions 

We present an improved, quasi-time-dependent radiative 
transfer model for the dust shell of the type IT A OH/IR 
star OH 104.9+2.4. This model is constrained by recent 
speckle-interferometric observations at 2.12/im and by 
spectro-photometric observations. The speckle observa- 
tions were carried out with the SAO 6 m telescope at 
a diffraction- limited resolution of 74 mas, which is suffi- 
cient for resolving the CDS at this wavelength. The re- 
constructed 2-dimensional visibility shows no major devi- 
ation from spherical symmetry. As a reference SED, we 
used the ISO SWS and LWS spectra (A = 2.38 - 196 fan) 
from 1996, which correspond to minimum pulsation phase, 
as well as additional photometric data at different phases 
(see Paper I for details). 

Our simultaneous modeling of visibilities at different 



pulsation phases and of the ISO SED of OH 104.9+2.4 
lead to modifications of our previous best-fitting model 
presented in Paper I. The new results demonstrate the 
need to consider changes in bolometric flux, as well as in 
optical depth and central star effective temperature with 
pulsation phase, in order to obtain more reliable results. 

The final best-fitting model of this study leads to pre- 
diction of a temperature at the inner rim of the dust shell 
of Tj n = 1000 K. The dust shell has a relative thickness of 
10 5 , while the density decreases outwards according to a 
radiatively driven wind as described in Ivezic et al. (fT999 ) . 
The dust grains follow a standard MRN distribution with 
grain sizes a = 0.005 — 0.28 /j,m, while the optical constants 
of the dust itself are represented best by the cold silicates 
from Suh <{T39"9"|l . The optical depth varies between min- 
imum and maximum pulsation phase from T2.2^m = 8.5 
to 7"2.2/Lim — 3.5, while the bolometric flux changes from 
Fboi = 0.7 x 10- 10 W/m 2 to F hol = 2.3 x lO" 10 W/m 2 . 

Based on pulsation models for Mira stars (Bessell et al. 
1996), variation of the effective temperature of the central 
star with pulsation phase is taken into account. For the 
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Fig. 5. Variation of the SED with pulsation phase. All 
data collected in Paper I (see Fig. 2 in Paper I and refer- 
ences therein) are plotted in order to test the quality of 
the SED fit of M-II for the 2002/2003 visibility phases. 

cool star OH 104.9+2.4, the best-fitting model predicts a 
change from T e g = 2250 K to T c g = 3150 K between min- 
imum and maximum pulsation phases. While a change in 
T e ff shows only a minor impact on SED and/or visibility, 
it leads to a prediction of the change in the stellar radius 
(from R+ = 729 Rq at minimum phase to i?* = 675 R@ 
at maximum phase) for our best-fitting model, which is in 
agreement with observations of similar objects (see, e.g., 
Bessell et al. 1995J ■ The inner rim of the dust shell moves 
outwards by roughly a factor of two from minimum to 
maximum phase (i?j n = 28.3 — 55.0 AU), while an increase 
of the mass-loss rate from M — 3.09 x 10~ 5 M Q /yr to 
Al = 5.69 x lO _5 M0/yr is found. Several observational 
constraints (1996 ISO SED, 2002/2003 SAO visibilities, 
IR photometry at different epochs) are represented very 
well by this model. Only the 2003 visibility at high spatial 
frequencies shows deviation from the model, which may 
be due to underestimation of observational errors. 

In the future, it would be very important to obtain K'- 
band visibility measurements at pulsation phase $ = 0.5 
in order to check the validity of the phase-dependent 
model, although conclusions may be limited due to pos- 
sible cycle-to-cycle variations. Additional measurements 
for intermediate phases would help to check if the opti- 
cal depth really scales linearly with the pulsation ampli- 
tude or if a more complex approach has to be considered. 
Likewise, it is still necessary to carry out more photo- 
metric measurements in order to pin down the pulsation 
period with higher precision. Furthermore, visibilities at 
other wavelengths (probably J or H, which poses a chal- 
lenge due to the detection limit) would help to improve 
the modeling and would give more detailed information 
on the spatial structure of the dust shell. 
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